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Abstract 

We present an in depth discussion of the production of gravitational waves from an inflationary phase 
that could have occurred in the early universe, giving derivations for the resulting spectrum and 
energy density. We also consider the large-scale anisotropy in the cosmic microwave background 
radiation coming from these waves. Assuming that the observed quadrupole anisotropy comes 
mostly from gravitational waves (consistent with the predictions of a flat spectrum of scalar density 
perturbations and the measured dipole anisotropy) we describe in detail how to derive a value for 
the scale of inflation of (1.5 — 5) x 10 16 GeV, which is at a particularly interesting scale for particle 
physics. This upper limit corresponds to a 95% confidence level upper limit on the scale of inflation 
assuming only that the quadrupole anisotropy from gravitational waves is not cancelled by another 
source. Direct detection of gravitational waves produced by inflation near this scale will have to 
wait for the next generation of detectors. 
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1. Introduction 

It has long been realized that a period of inflation in the early universe would lead to production 
of a well defined spectrum of gravitational waves |T[ (while inflation is not the only method of gen- 
erating a stochastic background of horizon-sized waves 0, it is certainly the most well motivated). 
Since the waves decouple early they are a potentially good probe of conditions in the early universe. 
In fact as pointed out in || |[ [7| such waves, which are fluctuations in the metric, can result 
in distortions of the cosmic microwave background radiation (CMBR) thus allowing information 
on the amplitude of the waves (and the parameters of the inflationary phase that produced them) 
to be derived from the measured CMBR anisotropy. This program proceeds in 3 steps. First the 
expected fluctuation spectrum is derived, then the consequences of such a spectrum for CMBR 
anisotropies are calculated. Finally comparing these predictions with the observations allows us to 
infer the parameters of an inflationary theory that produces such waves. 

Recently an analysis of the implications of a CMBR anisotropy connected with gravitational 
waves from inflation was presented ||. In this paper we fill in details of the results presented in 
II and further examine the issues addressed therein. A review of each of the 3 steps mentioned 
above will be given (many of the analytic results presented here have appeared in some form 
scattered in the literature, we have attempted to check, unify and reconcile previous estimates). 
The results can be compared to the recent COBE data || which gives the first positive measurement 



of anisotropies in the CMBR (previously this anisotropy had only been limited from above [|T^, 11 1). 
On the theoretical side we give the details of derivations of some of the key results in the hope 
that these will be of use. We present in one place the derivation of the scale invariant spectrum 
predicted by exponential inflation (and a description of the scalar analogy), the energy density 
p or equivalently Q g , the CMBR anisotropy produced and the predicted temperature correlation 
function. We show how "redshift during horizon crossing" dramatically reduces the predicted Q g 
from their asymptotic value for waves re-entering the horizon during the matter dominated epoch. 
Our results for the multipole moments of the CMBR temperature fluctuation agree with [Q (though 
we use an approach due to |12 to derive the scale invariant spectrum) and are a factor of ~ 2 larger 
than Q . With these results in hand we describe how the recent COBE observation is consistent with 
inflation and can limit the "scale near the end of inflation" assuming that the CMBR anisotropy is 
due, at least in part, to inflation-produced gravitational waves (the contribution from scalar density 
perturbations is constrained by the dipole anisotropy as we shall discuss). The scale turns out to 
be fewxl0 16 GeV |8| - an interesting energy for particle physics. 

The outline of the paper is as follows: section 2 establishes the notation and conventions. In 
sections 3 and 4 we discuss the scale-invariant spectrum predicted by exponential inflation and the 
spectral energy density. Section 5 is devoted to the derivation of the anisotropy generated by a 
stochastic spectrum of gravitational waves. The relation to the temperature correlation function 
measured by double- and triple-beam experiments (e.g. COBE) is derived in section 6 and the 
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current measurement of the CMB anisotropy is compared to the theoretical predictions in section 
7. Section 8 contains the conclusions. 



2. Cosmological history 

We will be interested in a cosmology which had an early period of inflation (where the dominant 
energy density resided in the vacuum). To be specific let us suppose that the universe underwent 
exponential inflation, then became radiation dominated and is currently in a matter dominated 
phase. The assumption of exponential inflation is a good approximation for "successful" inflation 
models based on general relativity since the requirement of a "slow roll" period in the evolution 



of the inflaton field requires a very flat potential. The other possibility is that of power law inflation, 
where the scale factor grows like a power of time. This changes the energy density spectrum (putting 
more power in long wavelength modes) and enhances the CM BR anisotropies for small I [f| . It seems 
difficult to have consistent power law inflation even in Brans-Dicke theory and we will not consider 
it further. 

Our metric is of the usual k = Robertson- Walker form 

ds 2 = -dt 2 + R 2 {t) (dx 2 + dy 2 + dz 2 ) = R 2 (t) (-dr 2 + dx 

where dr = dt/R(t) is the conformal time. For an equation of state p = qp the scale factor is a 
power of the conformal time, R(t) ~ r 2 ^ 1+3q \ For later convenience we list below the scale factor, 
which we normalize to unity today, for the history described above along with the dominant form 
of energy density for each epoch. We assume that the transitions are sudden and match R(r) and 
R{t) at the transition points (this will be a good enough approximation for our purposes) jJEJ. This 
matching implies that r is discontinuous across each transition and R(r) has the form 

{—(Ht)^ 1 vacuum rG (— oo, — r 2 ) 
It^It 2 radiation rG (r 2 ,r 1 /2) (2) 
t 2 /tq matter rG (ti,t ) 

where H is the Hubble constant during inflation. Equating /J* dr with J X dt/R(t) we find r = 3t 
is the conformal time today and T\ = (SGtitg) 1 ^ 3 is the conformal time at matter-radiation equality 
(both To and T\ refer to r as measured in the matter dominated era - in this respect our notation 
differs from that in [ 12 1 ) . The conformal time at the end of inflation is r 2 = r / yj 2Ht\ . The Hubble 



constant H and vacuum energy density Vq driving the inflation are related by 

2 _ 8tt Vq _ 8vr , 

H ~ " Y mpiV (3) 

where v = Vo/m^. 
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It is useful to have an expression for the size of a wave, relative to the horizon size. During a 
normal FRW expansion the coordinate radius of the horizon grows as dr = dt/R(t) = dr. If the 
period of inflation is long the coordinate radius of the horizon is very nearly zero at the beginning 
of the radiation dominated era. Thus at time i* the proper radius of the horizon is ~ i?(r*)r*. If we 
consider a wave just entering the horizon at this time we find its comoving wavenumber, k, must 
satisfy 

Aphys = R*T* => k = - = (4) 

Aphys T* 

This leads us to the following two limits. A wave with kr <C 2tt is well outside the horizon while a 
wave with kr ^> 2n is well within the horizon. 

3. The scalar analogy 

A classical gravitational wave in the linearized theory is a ripple on the background space-time 

9nv = R 2 {t) (r)^ + h^) where r)^ = diag(-l, 1, 1, 1), < 1 (5) 

In what follows we will work in transverse traceless (TT) gauge and denote the two independent 
polarization states of the wave as +, x. In the linear theory the TT metric fluctuations are gauge 
invariant (they can be related to components of the curvature tensor) ||15|| . We can write a plane 
wave with comoving wavenumber \k\ 

h^ u (r, x) = h x (r; k)e ik - 3 e^ u (k; A) (6) 

where e pv (k\ A) is the polarization tensor and A = +, x. The equation for the amplitude h\(r; k) is 
obtained by requiring the perturbed metric (0) satisfy Einstein's equations to 0(h). One finds [jnj 

h\ + 2^h x + k 2 h x = (7) 



As noticed by Grishchuk this is just the massless Klein-Gordon equation for a plane wave in the 
background space-time. Thus each polarization state of the wave behaves as a massless, minimally 
coupled, real scalar field, with a normalization factor of VTQttG relating the two. We can use the 
following heuristic argument to fix the y/lQirG: in the linearized theory one can write the Hilbert 
action for the fluctuation using the contribution to the Ricci scalar from h^. From the derivation 
of the stress-energy tensor T$ @ one sees = \h^ u ^ p h^ u ' p so the corresponding Hilbert action 
is 



(/,) 



S H 




(Vh + f + [Vhy 
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which is the action for two real, massless, scalar fields +jX = (lQixG)~ 1 l 2 h + ^ x as expected. Thus 
the study of quantum mechanical graviton production reduces to the study of the fluctuations of a 
scalar field in the curved background space-time, a factor \^16nG relating the two cases. 

4. The spectrum and Vt g 

Quantum fluctuations have an important consequence in a cosmology with inflation. During 
inflation short wavelength quantum fluctuations will get red-shifted out of the horizon, after which 
they freeze in. In terms of (0) the freezing in of extrahorizontal modes is the statement that when 
wavelength is much larger than the horizon the k 2 h term is negligible and the solution is a constant 
amplitude: h = (even before this limit is reached we note that the behaviour of waves outside the 
horizon is qualitatively different from the damped oscillatory behaviour of waves within it). When 
the mode re-enters the horizon at a much later epoch it appears as a long wavelength, classical 
gravitational wave (in analogy with the case of scalar fluctuations considered by [JT/j|). One way of 



thinking of this |18[ is that the number of quanta describing a state of constant amplitude grows 
as a power of the scale factor, the energy of each quantum being redshifted. Thus when the mode 
re-enters the horizon it represents a very large number of quanta. 

The spectrum of gravitational waves generated by quantum fluctuations during the inflationary 
period can be derived by a sequence of Bogoliubov transformations relating creation and annihilation 
operators defined in the various phases: inflationary, radiation and matter dominated |L2], [19| . 



The key idea is that for modes which have inflated outside the horizon the transitions between 
the phases are sudden and the universe will remain in the quantum state it occupied before the 
transition (treating each of the transitions as instantaneous is a good approximation for all but the 
highest frequency graviton modes). However the creation and annihilation operators that describe 
the particles in the state are related by a Bogoliubov transformation, so the quantum expectation 



value of any string of fields is changed (see e.g. |TJ| for a discussion of the calculation of Bogoliubov 
coefficients). 

We calculate the statistical average of the ensemble of classical waves by considering the corre- 
sponding quantum average. The simplest examples to consider are the 2-point functions. For the 



quantum theory we calculate the scalar 2-point function, following JT2 

k 3 

W) 



A = [ d 3 x e ils (ip\ <f)(x, r WO, r) U) where U) = Ide Sitter vac.) (9) 

(27rr J 



with the field 

(£, r ) = J d 3 k (e^a(A;)0 fc (r) + e^V^^r)) (10) 
where 4>k{ T ) is a properly normalized solution of (0). For the case of interest, a matter dominated 



5 



universe today, 



ifc is related to a Hankel function of order 3/2 

„— ikr 



2k (2tt) 3 / 2 R(t) 



kr 



In the exponentially inflating phase has the same functional form (with different R(r)) and in 
the radiation dominated phase the form is again as (0), with appropriate R{t) and also with the 
factor in parentheses absent. 

The quantum 2-point function is obtained by using the Bogoliubov coefficients relating the 
creation and annihilation operators a and a) of the field in the 3 phases 



^rad 
Q"mat 



ci(k)a inf (k) + c* 2 {k)a\ nf {-k) 
c 3 (k)a in f(k) + c\(k)a\ nf (-k) 



(12) 
(13) 



We are interested in waves which are still well outside the horizon at the time of matter-radiation 
equality (kr\ <C 2n) since these will give the largest contribution to the CMBR anisotropy today. 
Matching the field and its first derivative at t 2 , T\ in the limit kr <C 2n we find 



Ci « -c 2 



-Ht x 
(kr o y 



c 3 



-3iH 
2kHl 



(14) 



which is, correcting for differences in conventions, in agreement with [12]. Using (|Tl|Jl3| Jl 
one obtains, for waves re-entering the horizon in the matter dominated era 



in 



H 2 



2(2nf 



3ji(fcr) 
kr 



A 



H 2 



GW 



PI 



3ji(fcr) 
kr 
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(15) 



where in the last step the scalar result was multiplied by IQnG to obtain the corresponding result for 
gravitational waves. The corresponding expression for waves entering during the matter dominated 
era can be obtained by replacing the factor in brackets by jo(kr) as one would expect (since jo gives 
the right r dependence from (]?]) and the amplitude as kr — > should be the same). 
We match this to our classical ensemble of gravitational waves, h\{r; k), by writing 



h x {r; k) = A(k)a x {k) 



3ji(fcr) 
kr 



with A = +, x 



(16) 



where the term in [• ■ ■] is the real solution to (0) in the matter dominated phase and a\(k) is a 
random variable with statistical expectation value 



a x (k)ay(q)) = k-^ 3 \k~q)S X y 



(17) 
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The 2-point function analogous to A in this case is simply 



A', 



AA' 



A 2 {k) 



3ji(fcr 
kr 



1 2 
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AA' 



Matching the quantum and classical 2-point functions gives us the well known prediction for the 
(fc-independent) spectrum of gravitational waves generated by inflation 



A\k) 



H 2 



PI 



3vr 



(19) 



The reason for the explicit factor of k 3 introduced in the definition of A (|9]) can be seen by 
considering the expression for the spectral energy density in (classical) gravitational waves. In our 
notation 



P 



k 2 

phys 

16nG 



tii + h 2 



(20) 



where (• • •) indicates an average over many wavelengths/periods as well as an average over the 
stochastic variable a\{k). k p h ys = k/R(r) is the physical wavenumber. For a fixed (non-stochastic) 
spectrum of waves we could write the average over space as an integral over Fourier modes using 
Parseval's theorem and then relate the "power spectrum" \h +jX (k)\ 2 to the correlation function (§). 
With (|T6|) and (|17D as defined the stochastic average produces the same result. Explicitly 



k 



dpx 
dk 



U2 

n phys "XT 

AG 



A = +, x 



(21) 



which goes as R 4 (r) as expected (an overbar denotes averaging). Up to the time evolution factor 



(3ji(kr)/kT) 2 for matter and jo(kr) for radiation dominated phases this agrees with [pO |. Including 
this factor for waves re-entering during the matter dominated phase (i.e. considering evolution as 
the wave enters the horizon) leads to a significant suppression of p. Writing the energy density as 
a fraction of the closure density one finds, for waves just entering the horizon M, 



a, 



yy k_ 

A=+,x Pc 



dk 



k=2ir/T 




RD 
MD 



(22) 



In deriving this use has been made of the fact that k p h ys = 1(2)ttHhc for waves crossing the horizon 
during the matter (radiation) dominated phase (Hhc is the Hubble constant at the time the wave 
re-enters the horizon). 

Will we be able to see this energy density with terrestrial or astrophysical gravitational wave 
detector? The most sensitive gravitational wave limit currently comes from the timing of millisecond 
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pulsars. This is sensitive to "short" wavelengths (periods of order years) and puts an upper limit 
on the energy density in such waves Q g < 9 x 10~ 8 (68%) (this limit improves as the fourth power 
of the observing time: T^ bs ) PB| . Such short wavelength modes entered the horizon during the 
radiation dominated era. If we write the Hubble constant today H = 100/i km/s/Mpc then these 
waves contribute to fl g today an amount suppressed by a factor of p ra d/ Pc ~ 4 x 10~ 5 /i~ 2 compared 
to horizon crossing since they redshift with one extra power of R compared to matter. For values of 
v ~ 10~ 9 , consistent with the anisotropy of the CMBR (see later), the energy density per logarithmic 
frequency interval in short wavelength modes is predicted to be Q g ~ 10~ 13 ||, which is ~ 6 orders 
of magnitude lower than the millisecond pulsar limit. Even the proposed LIGO detector |2{| which 



has a sensitivity of Q g ~ 10~ n is short of the mark by 2 orders of magnitude. Thus we expect that 
a positive detection of the signal will require a significant advance in technology. 

Since the predicted CMBR anisotropy multipoles, up to I ~ 9, from gravitational waves and 
scale invariant scalar fluctuations are very similar and direct detection of gravitational waves 
is still some way off the outlook for determining unambiguously that a significant fraction of the 
observed CMBR anisotropy comes from gravitational waves seems bleak. 

5. Anisotropy from the stochastic background 

It is conventional to expand the CMBR temperature anisotropy in spherical harmonics 

^(M) = E a '-MM) (23) 

Ira 

We can calculate the prediction of a given spectrum of gravitational waves in terms of the a/ m . The 
temperature fluctuation due to a gravitational wave can be found in the linearized theory to be 



[T23 

f =4I dA ^#^ (24) 

where A is a parameter along the unperturbed path and the lower (upper) limit of integration 
represents the point of emission (reception) of the photon. 

Now we project out a multipole and calculate the rotationally symmetric quantity 

ST ,„JT 




E / *Mif Y^YUn') /— (n) d T (n')\ (25) 



where using (|,|17].^4|) we have for a spectrum of waves 
/8T,„JT,„A 1 ^ rd 3 k 



dM /i(r';A;)e^'e^(fc;A)x' M x' 



(26) 
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and an overdot represents a derivative with respect to r. We use the rotational symmetry of af to 
make the replacement 

£ Y* m {n)Y lm {n') = £ Y lm(ttkx)Ylm(tt' kx ) (27) 



where Q kx indicates that the angles are defined with respect to k. We can evaluate directly 



(k; \)x^x v = (5f cos(20) + 5* sin(20)) sin 2 9 



where 8, <fi are the usual spherical angles relating k and x. 

First concentrate on the angular integrations. Expanding the exponential 



ik-x 



J2i n (2n + l)j n (kx)P n (cos9, 



kx ) 



n=0 



it is not hard to show (formula 7.125 of Gradshteyn and Ryzhik |M is useful) that 



hm(k, x) 



7T, 



dn kx Y; m (n kx ) ^^(k- 

21 + 1(1 + 2)! 
\ 4tt (Z-2)! 



H £ i n (2n + l)j n (kx) (c_ 2 ^_ 2 + CoK + c 2 5? +2 



where 



l/C-2 : 

1/co = 

l/c +2 = (2/ + l)(2/ + 3)(/ + 5/2) 
Inserting this back into the expression for af and using (JT 



in 



(^ + -^A X )e 2 + (^+^A X H 

(21- l)(2Z + l)(Z-3/2) 
--(2Z - 1)(2Z + 3)(Z + 1/2) 



E E 

m A=+,x 



dkdQk 
k 



dA A(k) 



d /3ji(Zcr) 



cZr 



Ilm(k, X / 



(21 



(29) 



(30) 



(31) 
(32) 

(33) 
(34) 



(35) 



The integral over the path dA can be parameterized by the distance from the origin along the 
line of sight, r, so \x(r)\ = r and r(r) = r — r. We will defer consideration of the limits of the k 
integral for the moment. Thus 



af) = 3Qn 2 (2l + l 



(j+2)! 
(Z-2)! 



kdk A 2 (k)\F l (k)[ 



(36) 
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where the function Fi(k) is defined as 

Ji-2(kr) 2j l (kr) Jwikr) 



Fi(k) 




(2Z-1)(2Z + 1) (2Z — 1)(2Z + 3) (2Z + l)(2Z + 3) 



(37) 



Accounting for the factor of two difference between the definitions of A 2 (k) this is precisely the 
result of ||, and is a factor of « 2 larger than the earlier result of Q. The integral in equation 
(3"7p can be expressed in terms of elementary functions and the sine integral but the result is very 



cumbersome and is not presented here. 

For waves entering the horizon during the matter dominated regime the upper limit of the k 
integration is 2tt/ti. The results are very insensitive to the exact choice of upper limit since the 
integral is dominated by small k « 2-k/tq, i.e. waves that have recently entered the horizon. For 
very high frequency modes the transition between phases of the universe is no longer sudden and 
graviton production is suppressed. Taking this high frequency cutoff to be k = 2tt/ti (i.e. restricting 
attention to waves that entered the horizon during the matter dominated era) introduces negligible 
error. The lower limit of the k integral is zero, since waves of arbitrarily long wavelength can 
contribute^. From the form of the integrand however one can see that the contribution to (af) 
tends to zero as k n (n > 1) as kr — > 0. Thus the contribution from very long wavelength modes is 
suppressed as one might expect on physical grounds. 

One can now calculate the predicted (af) for any spectrum A(k) of gravitational waves from 
(PB|). For the scale-invariant spectrum ( |I9"D of exponential inflation the (af) are shown in table 
The calculation of the expectation value (af) is not the end of the story however. Before we compare 
predictions with observations we must also consider the statistical properties of af. The fact that 
the "weakly coupled" nature of inflation predicts the a\ m to be independent makes this problem 
tractable. Given that each of the a\ m are independent Gaussian random variables the probability 
distribution for each af, whose mean (af) we calculated, is of the \ 2 form and can be written 

i-i/2 e -y 21 + 1 af 

p (y)dy = —t, — r\ dy where v = o / 2\ ( 38 ) 

I t + o Z \ a l) 



2) 

This agrees with |J. One can calculate the confidence levels for af in terms of the incomplete 
gamma function (values for the 68, 90 and 95% (lower) confidence levels can be found in table [l]). 
We note in passing that the modal value of af for any universe is {21 — 1)/(2Z + 1) times the mean. 

6. The correlation function 

Limits quoted on the multipole moments af can be directly compared to (|36|) to constrain the 
inflationary theory. However many experiments (and importantly COBE) are of the double- or 

2 We thank Mark Wise and Vince Moncrief for very helpful discussions on this point. 
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triple-beam type and do not simply measure any one multipole moment but can constrain the 
correlation function || |20 



C(9 



'5T 



5T 



21', a J 



(39) 



21 



where the average is over all positions x\, x 2 on the sky with x\ ■ x 2 fixed, and a is a measure of the 
angular response of the detector. We shall take the angular response to be Gaussian [|21] 

e 2 



dR(6,~ 



exp 



2a 2 



(40) 



where 8, are the angles relative to the beam direction and o <C 1 which justifies the use of the 
small angle expansion. We define the "smeared" ST/T by 



T 



[xi, a 



5>, m | dQ' R(Q')Y lm (Q) 

lm 

lm ^° 



dx Pi(l - a 2 x)e x )Y lm (x l ) 



(41) 
(42) 



where Q' is the direction Q measured relative to x\ and we have approximated cos#' « 1 — 6' 2 /2 
for small 9'. Note that the term in parenthesis suppresses the higher I modes as we would expect 
physically. We can expand this term for small a 2 



dx Pi{l-a 2 x)e~ x = ^(-a 2 ) n P/ n) (l) 

?1=0 



n=0 



—a 



J_(H~n)! 
' 2 n n\ (I - n)\ 



(43) 
(44) 



Making the approximation (I + n)\ / (I — n)\ « (/ + l/2) 2n for large I we can write the suppression 
term as: exp [— (/ + l/2) 2 o" 2 /2] (and at the same level of approximation we can replace (I + 1/2) 2 
with 1(1 + 1) as used by This exponential suppression proves a good approximation for o < 0.1 
which is the range for most detectors. A simple exercise in manipulating spherical tensors allows 
us to write 



{Y lm {x 1 )Y l l m> {x 2 )), 



1 



5lV $mm> Pl(xi ■ X 2 ) 



' 21 4vr 

so we can express the correlation function ([39|) in a simple from as (see || |20[| ' 

1 oo 

C ( Xl . X2 - a) = ± £ a 2 P l ( Xl ■ x 2 ) e^ l+1 l 2 ^ 2 
47r 1=0 



(45) 



(46) 
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The predicted correlation function using the moments of table [I] is shown in figure |1] for a 
gaussian beam width of o = 10°. The error bars represent the upper and lower limit of the 68% 
confidence region arrived at using for the distribution of af/v (there is no simple analytic form 
for the probability distribution of C(6,a) so figure [I] was obtained from a Monte-Carlo) and the 
solid line is C(9,a) evaluated using af = (af). The correlation function inferred from the COBE 
sky maps || is in agreement with the form shown in figure |l|. 

7. Comparison with observations 

The recent results of COBE can be summarized as a non-zero value of the quadrupole moment af . 
If we require a Harrison-Zel'dovich spectrum the quadrupole, from the Q rm s-ps value, is measured 
to be |, 1 

a\ = (4.7 ±2) x 1(T 10 (47) 

An early inflationary phase produces not only gravitational waves but also a scale invariant 
spectrum of scalar density perturbations. To limit the contribution of these latter to a\ we can 
require the induced dipole from long wavelength scalar modes not greatly exceed the observed 
dipole anisotropy. For a given (e.g. flat) spectrum one can use this to place limits on all the 
higher multipoles. At the 90% confidence level an upper limit of a\ ~ 2 x 10~ 10 has been 
derived (ignoring the transfer function corrections). Equivalently fitting the observed clustering to 
a primordial spectrum yields a 2 , ~ (2 — 10) x 10~ n p5f . Such estimates suggest a large fraction of 
the quadrupole anisotropy may be due to long wavelength gravitational waves. 

Let us assume that the observed quadrupole is due entirely to gravitational waves. Using {a^/v 
from table [l] and including the distribution for the measured a\ (gaussian) and predicted a\/v (|38|) 



one can infer a distribution for v. The result obtained from a Monte-Carlo is shown in figure 0. 
This corresponds to a mean v — 6.1 x 10~ n , and a modal v = 4 x 10~ n with non-gaussian errors. 
The allowed range of v from figure ^] can be summarized as [H] 

2.3 x 10- 11 < v < 1.5 x 10~ 10 68% CL . . 

2.5 x 10~ 12 < v < 3.7 x 10- 10 95% CL ^ > 

One can scale these results directly to obtain values for v if gravitational wave contribution to a\ is 
not 100%, the rest being made up or shielded by scalar fluctuations for example. 

Converting to a more familiar energy scale M = (f 1//4 Mp^) the 95% confidence level range 
above is M E (1.5 - 5.2) x 10 16 GeV, with a best fit M = 3 x 10 16 GeV. The upper limit quoted 
corresponds to a strict upper limit on the scale of inflation assuming that the quadrupole moment 
due to gravitational waves is not being cancelled to any significant degree by other sources (e.g. 
scalar density perturbations). A large cancellation between two sources of anisotropy would be 
unlikely. 
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No realistic model of inflation currently exists, however most are based on the idea that a phase 
transition occurred in the early universe to produce the vacuum energy density to drive the inflation. 
If such a transition were to occur then the natural scale would seem to be the scale of unification 
of the coupling constants in the context of a Grand Unified Theory (GUT) (an exception is the 
case of chaotic inflation). The recent precision measurements from LEP have allowed us to sharpen 
our predictions for this scale. Assuming a supersymmetric GUT with the supersymmetry breaking 
scale ~ ITeV, one predicts [26] a unification scale of (1 — 3.6) x 10 16 GeV. For a theory based 
on SO(10) one finds M = 10 15 ' 8± ' 22 GeV [27]. The limits from proton decay experiments are very 
theory dependent, however they give a lower bound on M. The current limits are near or larger 
than 10 15 GeV for fashionable GUT theories. 

Unless there is significant fine tuning or hierarchies in the theory the vacuum energy density 
associated with a transition from a GUT with scale Mq is Vq = (kMg) 4 where k is a (calculable) 
number usually in the range k 4 = .01 — 1. This puts the GUT scale very close to the scale of 
inflation inferred above - a "numerical coincidence" which is both suggestive and exciting. 

8. Conclusions 

As the observations of CMBR anisotropies are refined and other gravitational wave detectors 
come on line the results presented here can allow a comparison of theoretical predictions and 
experimental results, which can help to constrain models of inflation and early universe physics. 

The recent COBE results allow the possibility that the observed quadrupole anisotropy results 
from gravitational waves generated during inflation with a scale of (1.5 — 5.2) x 10 16 GeV - a very 
interesting range of energies. If the anisotropy is due to gravitational waves, confirmation will 
probably not come from existing or currently planned gravitational wave detectors. Further results 
are eagerly awaited. 
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Figure Captions 



Figure 1: Results of a Monte-Carlo for the predicted correlation function C(6;a) generated by a 
flat spectrum of gravitational waves. The error bars represent the upper and lower limits of the 
68% confidence region arrived at using the predicted distribution of af/v. The solid curve is the 
correlation function calculated using (af). The beam width is taken to be a = 10°. 



Figure 2: The probability distribution for the scale v = Vo/m^ of inflation as determined by Monte- 
Carlo from the COBE measurement of the quadrupole anistropy assuming that the whole anistropy 
is due to gravitational waves. 



Tables 
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<# /(a?) 


/ 


(tf)/v 


68/90/95% CL 


2 


7.74 


.63/.32/.23 


3 


4.25 


.69/.40/.31 


4 


3.10 


J3/.46/.37 


5 


2.50 


.76/.51/.42 


6 


2.12 


J8/.54/.45 


7 


1.85 


.80/.57/.48 


8 


1.64 


.81/.59/.51 


9 


1.48 


.82/.61/.53 


10 


1.35 


.83/.63/.55 



Table 1: Multipole coefficients a; 2 for modes Z = 2 — 10 predicted for a stochastic background of 
gravitational waves generated by exponential inflation. The 68, 90 and 95% (lower) confidence 
levels are also shown as fractions of (af). 
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